{smcl}
{com}{sf}{ul off}{txt}{.-}
      name:  {res}<unnamed>
       {txt}log:  {res}/N/project/suicide_study/pnas_replication/results/log/mi_3.smcl
  {txt}log type:  {res}smcl
 {txt}opened on:  {res}20 Aug 2020, 05:05:11
{txt}
{com}. 
. if ("`model'" == "logit"){c -(}
.         use "${c -(}home_dir{c )-}/data/processed/suicide_reg_v1_raw.dta", clear
. {c )-}
{txt}
{com}. 
. if ("`model'" == "mi") {c -(}
.         use "${c -(}home_dir{c )-}/data/processed/suicide_reg_v1_imputed_M10.dta", clear  
. {c )-}
{txt}
{com}. 
. * for now, we use the following simple survey weights 
. if ("`model'" == "mi") {c -(}
.         mi svyset `geo_type' [pw=ObsWgt0] 
{res}
      {txt}pweight:{col 16}{res}ObsWgt0
          {txt}VCE:{col 16}{res}linearized
  {txt}Single unit:{col 16}{res}missing
     {txt}Strata 1:{col 16}<one>
         SU 1:{col 16}{res}county
        {txt}FPC 1:{col 16}<zero>
{p2colreset}{...}
{res}{com}. {c )-}
{txt}
{com}. else {c -(}
.         svyset `geo_type' [pw=ObsWgt0]  
. {c )-}
{txt}
{com}. 
. * margins for each category
. program margin_interact 
{txt}  1{com}.         args X Y k model
{txt}  2{com}.         sum `X', d 
{txt}  3{com}.         local gap = (`r(max)' - `r(min)') / `k' 
{txt}  4{com}.         if ("`model'" == "logit") {c -(}
{txt}  5{com}.                 margin `Y', at(`X' = (`r(min)' (`gap') `r(max)')) predict(pr)
{txt}  6{com}.         {c )-} 
{txt}  7{com}.         else if ("`model'" == "mi") {c -(}
{txt}  8{com}.                 mimrgns `Y', at(`X' = (`r(min)' (`gap') `r(max)')) predict(pr)
{txt}  9{com}.         {c )-}       
{txt} 10{com}. end 
{txt}
{com}. 
. program mchange_mi
{txt}  1{com}.         args X k model
{txt}  2{com}.         if ("`model'" == "logit") {c -(}
{txt}  3{com}. 
.                 if ("`k'" == "continuous") {c -(}
{txt}  4{com}.                         sum `X' if e(sample), d 
{txt}  5{com}.                         margin, at(`X' = (`r(min)' `r(max)')) post predict(pr)
{txt}  6{com}.                         mlincom  2 - 1, decimal(7) stat(all)    
{txt}  7{com}.                 {c )-} 
{txt}  8{com}.                 else if ("`k'" == "binary") {c -(}
{txt}  9{com}.                         sum `X' if e(sample), d 
{txt} 10{com}.                         margin, at(`X' = (`r(min)' `r(max)')) post predict(pr)
{txt} 11{com}.                         mlincom  2 - 1, decimal(7) stat(all)
{txt} 12{com}.                 {c )-} 
{txt} 13{com}.                 else if ("`k'" == "categorical") {c -(}
{txt} 14{com}.                         margin `X' if e(sample)==1 , at() pwcompare predict(pr) 
{txt} 15{com}.                 {c )-}
{txt} 16{com}.         {c )-}
{txt} 17{com}. 
.         else if ("`model'" == "mi") {c -(}
{txt} 18{com}. 
.                 if ("`k'" == "continuous") {c -(}
{txt} 19{com}.                         sum `X' , d 
{txt} 20{com}.                         mimrgns, at(`X' = (`r(min)' `r(max)')) post predict(pr)
{txt} 21{com}.                         mlincom  2 - 1, decimal(7) stat(all)
{txt} 22{com}.                 {c )-} 
{txt} 23{com}.                 else if ("`k'" == "binary") {c -(}
{txt} 24{com}.                         sum `X' , d 
{txt} 25{com}.                         mimrgns, at(`X' = (`r(min)' `r(max)')) post predict(pr)
{txt} 26{com}.                         mlincom  2 - 1, decimal(7) stat(all)
{txt} 27{com}.                 {c )-} 
{txt} 28{com}.                 else if ("`k'" == "categorical") {c -(}
{txt} 29{com}.                         mimrgns `X' , at()   pwcompare predict(pr)      
{txt} 30{com}.                 {c )-}
{txt} 31{com}.         {c )-}
{txt} 32{com}. end
{txt}
{com}. 
. * create some dummy codings
. tab Race5, gen(race5_nh)

      {txt}Race5 {c |}      Freq.     Percent        Cum.
{hline 12}{c +}{hline 35}
          1 {c |}{res}107,831,845       79.53       79.53
{txt}          2 {c |}{res} 14,631,144       10.79       90.32
{txt}          3 {c |}{res}  2,867,304        2.11       92.44
{txt}          4 {c |}{res}  2,753,311        2.03       94.47
{txt}          5 {c |}{res}  7,499,129        5.53      100.00
{txt}{hline 12}{c +}{hline 35}
      Total {c |}{res}135,582,733      100.00
{txt}
{com}.         rename race5_nh1 White_nh 
{res}{txt}
{com}.         rename race5_nh2 Black_nh 
{res}{txt}
{com}.         rename race5_nh3 AIAN_nh 
{res}{txt}
{com}.         rename race5_nh4 AsPI_nh 
{res}{txt}
{com}.         rename race5_nh5 Hispanic
{res}{txt}
{com}. 
. tab MarStat5, gen(ms)

   {txt}MarStat5 {c |}      Freq.     Percent        Cum.
{hline 12}{c +}{hline 35}
          1 {c |}{res} 75,689,936       55.83       55.83
{txt}          2 {c |}{res} 10,166,664        7.50       63.32
{txt}          3 {c |}{res} 13,925,509       10.27       73.60
{txt}          4 {c |}{res}  2,609,308        1.92       75.52
{txt}          5 {c |}{res} 33,190,855       24.48      100.00
{txt}{hline 12}{c +}{hline 35}
      Total {c |}{res}135,582,272      100.00
{txt}
{com}.         rename ms1 Marrd5
{res}{txt}
{com}.         rename ms2 Widow5
{res}{txt}
{com}.         rename ms3 Divor5
{res}{txt}
{com}.         rename ms4 Separ5
{res}{txt}
{com}.         rename ms5 NvMar5
{res}{txt}
{com}. 
. tab AgeGrp4, gen(ag) 

    {txt}AgeGrp4 {c |}      Freq.     Percent        Cum.
{hline 12}{c +}{hline 35}
          1 {c |}{res} 21,231,562       15.66       15.66
{txt}          2 {c |}{res} 39,043,312       28.80       44.46
{txt}          3 {c |}{res} 47,351,260       34.92       79.38
{txt}          4 {c |}{res} 27,956,599       20.62      100.00
{txt}{hline 12}{c +}{hline 35}
      Total {c |}{res}135,582,733      100.00
{txt}
{com}.         rename ag1 Age_15_24
{res}{txt}
{com}.         rename ag2 Age_25_44
{res}{txt}
{com}.         rename ag3 Age_45_64
{res}{txt}
{com}.         rename ag4 Age_65_Up
{res}{txt}
{com}. 
. destring St, replace 
{txt}St: all characters numeric; {res}replaced {txt}as {res}byte
{txt}
{com}. 
. * set-up equations
. local religion Rat_GC_ProE Rat_GC_ProM Rat_GC_ProB Rat_GC_Cath Rat_GC_Jew Rat_GC_Oth
{txt}
{com}. local contextual_control Rat_Poverty Rat_Mig_Cum Pop_Den
{txt}
{com}. 
. local religion Rat_GC_ProE Rat_GC_ProM Rat_GC_ProB Rat_GC_Jew Rat_GC_Oth
{txt}
{com}. local contextual_control Rat_Poverty Rat_Mig_Cum Pop_Den
{txt}
{com}. 
. local demographics_raw i.Female c.RAT_Female i.AgeGrp4 c.RAT_AgeGrp4_2 c.RAT_AgeGrp4_3 c.RAT_AgeGrp4_4 i.Race5 c.RAT_Race5_2 c.RAT_Race5_3 c.RAT_Race5_4 c.RAT_Race5_5 i.BornUSA c.RAT_BornUSA i.MarStat5 c.RAT_MarStat5_2 c.RAT_MarStat5_3 c.RAT_MarStat5_4 c.RAT_MarStat5_5
{txt}
{com}. local demographics_same i.Female c.std_same_prop_Sex i.AgeGrp4 c.std_same_prop_AgeGrp4 i.Race5 c.std_same_prop_Race5 i.BornUSA c.std_same_prop_BornUSA i.MarStat5 c.std_same_prop_MarStat5 
{txt}
{com}. local demographics_inter i.Female##c.std_same_prop_Sex i.AgeGrp4##c.std_same_prop_AgeGrp4 i.Race5##c.std_same_prop_Race5 i.BornUSA##c.std_same_prop_BornUSA i.MarStat5##c.std_same_prop_MarStat5
{txt}
{com}. 
. 
. if (`model_version' == 1){c -(}
.         local model_eq i.Year `demographics_raw' `contextual_control' `religion' 
.         local margin_demographics Female RAT_Female AgeGrp4 RAT_AgeGrp4_2 RAT_AgeGrp4_3 RAT_AgeGrp4_4 Race5 RAT_Race5_2 RAT_Race5_3 RAT_Race5_4 RAT_Race5_5 BornUSA RAT_BornUSA MarStat5 RAT_MarStat5_2 RAT_MarStat5_3 RAT_MarStat5_4 RAT_MarStat5_5 
. {c )-}
{txt}
{com}. if (`model_version' == 2){c -(}
.         local model_eq i.Year `demographics_raw' `contextual_control' `religion' i.UnEmpl c.RAT_UnEmpl i.PhysProb c.RAT_PhysProb
.         local margin_demographics Female RAT_Female AgeGrp4 RAT_AgeGrp4_2 RAT_AgeGrp4_3 RAT_AgeGrp4_4 Race5 RAT_Race5_2 RAT_Race5_3 RAT_Race5_4 RAT_Race5_5 BornUSA RAT_BornUSA MarStat5 RAT_MarStat5_2 RAT_MarStat5_3 RAT_MarStat5_4 RAT_MarStat5_5 UnEmpl RAT_UnEmpl PhysProb RAT_PhysProb
. {c )-}
{txt}
{com}. if (`model_version' == 3){c -(}
.         local model_eq i.Year `demographics_same' `contextual_control' `religion' 
.         local margin_demographics Female std_same_prop_Sex AgeGrp4 std_same_prop_AgeGrp4 Race5 std_same_prop_Race5 BornUSA std_same_prop_BornUSA MarStat5 std_same_prop_MarStat5 
. {c )-}
{txt}
{com}. if (`model_version' == 4){c -(}
.         local model_eq i.Year `demographics_same' `contextual_control' `religion' i.UnEmpl c.std_same_prop_UnEmpl i.PhysProb c.std_same_prop_PhysProb
.         local margin_demographics Female std_same_prop_Sex AgeGrp4 std_same_prop_AgeGrp4 Race5 std_same_prop_Race5 BornUSA std_same_prop_BornUSA MarStat5 std_same_prop_MarStat5 UnEmpl std_same_prop_UnEmpl PhysProb std_same_prop_PhysProb
. {c )-}
{txt}
{com}. if (`model_version' == 5){c -(}
.         local model_eq i.Year `demographics_inter' `contextual_control' `religion' 
. {c )-}
{txt}
{com}. if (`model_version' == 6){c -(}
.         local model_eq i.Year `demographics_inter' `contextual_control' `religion' i.UnEmpl##c.std_same_prop_UnEmpl i.PhysProb##c.std_same_prop_PhysProb
. {c )-}       
{txt}
{com}. * test how long it would take.
. * mi estimate: svy: mean Suic 
. * local demographics i.Female c.RAT_Female i.AgeGrp4 c.RAT_AgeGrp4_2 c.RAT_AgeGrp4_3 c.RAT_AgeGrp4_4 i.Race5 c.RAT_Race5_2 c.RAT_Race5_3 c.RAT_Race5_4 c.RAT_Race5_5 i.BornUSA c.RAT_BornUSA i.MarStat5 c.RAT_MarStat5_2 c.RAT_MarStat5_3 c.RAT_MarStat5_4 c.RAT_MarStat5_5
. * mi estimate: svy: logit Suic i.St `demographics' UnEmpl RAT_UnEmpl PhysProb RAT_PhysProb
. 
. * main effects : margins
. if ("`model'" == "mi"){c -(}
. 
.         mi estimate: svy: logit Suic i.St `model_eq', or 
{res}
{txt}Multiple-imputation estimates{col 47}Imputations{col 65}= {res}          10
{txt}Survey: Logistic regression{col 47}Number of obs{col 65}= {res}  11,814,307

{txt}{col 1}Number of strata{col 19}= {res}        1{txt}{col 47}Population size{col 65}={res}  418,786,443
{txt}{col 1}Number of PSUs{col 19}= {res}      918
{txt}{col 47}Average RVI{col 65}= {res}      0.0006
{txt}{col 47}Largest FMI{col 65}= {res}      0.0046
{txt}{col 47}Complete DF{col 65}= {res}         917
{txt}DF adjustment:{ralign 15: {res:Small sample}}{col 47}DF:     min{col 65}= {res}      908.90
{txt}{col 47}        avg{col 65}= {res}      914.43
{txt}{col 47}        max{col 65}= {res}      915.01
{txt}Model F test:{ralign 16: {res:Equal FMI}}{col 47}F({res}  47{txt},{res}  915.0{txt}){col 65}= {res}      741.71
{txt}Within VCE type: {ralign 12:{res:Linearized}}{col 47}Prob > F{col 65}= {res}      0.0000

{txt}{hline 13}{c TT}{hline 11}{hline 11}{hline 9}{hline 8}{hline 13}{hline 12}
{col 1}        Suic{col 14}{c |}      Coef.{col 26}   Std. Err.{col 38}      t{col 46}   P>|t|{col 54}     [95% Con{col 67}f. Interval]
{hline 13}{c +}{hline 11}{hline 11}{hline 9}{hline 8}{hline 13}{hline 12}
{space 10}St {c |}
{space 10}8  {c |}{col 14}{res}{space 2}-.1214146{col 26}{space 2} .0992626{col 37}{space 1}   -1.22{col 46}{space 3}0.222{col 54}{space 4}-.3162235{col 67}{space 3} .0733943
{txt}{space 9}13  {c |}{col 14}{res}{space 2}-.4962537{col 26}{space 2} .0984382{col 37}{space 1}   -5.04{col 46}{space 3}0.000{col 54}{space 4}-.6894446{col 67}{space 3}-.3030628
{txt}{space 9}21  {c |}{col 14}{res}{space 2}-.6283719{col 26}{space 2} .1019798{col 37}{space 1}   -6.16{col 46}{space 3}0.000{col 54}{space 4}-.8285135{col 67}{space 3}-.4282304
{txt}{space 9}24  {c |}{col 14}{res}{space 2}-.6786501{col 26}{space 2} .0994859{col 37}{space 1}   -6.82{col 46}{space 3}0.000{col 54}{space 4}-.8738972{col 67}{space 3}-.4834031
{txt}{space 9}25  {c |}{col 14}{res}{space 2}-.9494672{col 26}{space 2} .0998654{col 37}{space 1}   -9.51{col 46}{space 3}0.000{col 54}{space 4}-1.145459{col 67}{space 3}-.7534753
{txt}{space 9}34  {c |}{col 14}{res}{space 2}-.8417754{col 26}{space 2} .1008358{col 37}{space 1}   -8.35{col 46}{space 3}0.000{col 54}{space 4}-1.039672{col 67}{space 3}-.6438791
{txt}{space 9}35  {c |}{col 14}{res}{space 2} .0113567{col 26}{space 2}  .101985{col 37}{space 1}    0.11{col 46}{space 3}0.911{col 54}{space 4} -.188795{col 67}{space 3} .2115084
{txt}{space 9}37  {c |}{col 14}{res}{space 2}-.4502756{col 26}{space 2} .1024941{col 37}{space 1}   -4.39{col 46}{space 3}0.000{col 54}{space 4}-.6514263{col 67}{space 3}-.2491248
{txt}{space 9}40  {c |}{col 14}{res}{space 2}-.6219542{col 26}{space 2} .1075826{col 37}{space 1}   -5.78{col 46}{space 3}0.000{col 54}{space 4}-.8330914{col 67}{space 3}-.4108169
{txt}{space 9}41  {c |}{col 14}{res}{space 2}-.2975567{col 26}{space 2} .0997339{col 37}{space 1}   -2.98{col 46}{space 3}0.003{col 54}{space 4}-.4932904{col 67}{space 3}-.1018229
{txt}{space 9}44  {c |}{col 14}{res}{space 2}-.8165612{col 26}{space 2}  .099305{col 37}{space 1}   -8.22{col 46}{space 3}0.000{col 54}{space 4}-1.011453{col 67}{space 3}-.6216692
{txt}{space 9}45  {c |}{col 14}{res}{space 2}-.4875654{col 26}{space 2} .0994751{col 37}{space 1}   -4.90{col 46}{space 3}0.000{col 54}{space 4}-.6827913{col 67}{space 3}-.2923395
{txt}{space 9}49  {c |}{col 14}{res}{space 2}  .566737{col 26}{space 2} .1784937{col 37}{space 1}    3.18{col 46}{space 3}0.002{col 54}{space 4} .2164324{col 67}{space 3} .9170417
{txt}{space 9}51  {c |}{col 14}{res}{space 2}-.3606947{col 26}{space 2} .0967352{col 37}{space 1}   -3.73{col 46}{space 3}0.000{col 54}{space 4}-.5505434{col 67}{space 3} -.170846
{txt}{space 9}55  {c |}{col 14}{res}{space 2}-.6354016{col 26}{space 2} .1015165{col 37}{space 1}   -6.26{col 46}{space 3}0.000{col 54}{space 4}-.8346337{col 67}{space 3}-.4361695
{txt}{space 12} {c |}
{space 8}Year {c |}
{space 7}2006  {c |}{col 14}{res}{space 2}-.0618303{col 26}{space 2} .0181075{col 37}{space 1}   -3.41{col 46}{space 3}0.001{col 54}{space 4}-.0973675{col 67}{space 3}-.0262932
{txt}{space 7}2007  {c |}{col 14}{res}{space 2} .0145655{col 26}{space 2} .0183271{col 37}{space 1}    0.79{col 46}{space 3}0.427{col 54}{space 4}-.0214026{col 67}{space 3} .0505336
{txt}{space 7}2008  {c |}{col 14}{res}{space 2} .0160904{col 26}{space 2} .0200133{col 37}{space 1}    0.80{col 46}{space 3}0.422{col 54}{space 4} -.023187{col 67}{space 3} .0553677
{txt}{space 7}2009  {c |}{col 14}{res}{space 2} .0615614{col 26}{space 2} .0187613{col 37}{space 1}    3.28{col 46}{space 3}0.001{col 54}{space 4} .0247411{col 67}{space 3} .0983816
{txt}{space 7}2010  {c |}{col 14}{res}{space 2} .0446807{col 26}{space 2} .0203698{col 37}{space 1}    2.19{col 46}{space 3}0.029{col 54}{space 4} .0047038{col 67}{space 3} .0846576
{txt}{space 7}2011  {c |}{col 14}{res}{space 2} .0789179{col 26}{space 2} .0220904{col 37}{space 1}    3.57{col 46}{space 3}0.000{col 54}{space 4} .0355642{col 67}{space 3} .1222717
{txt}{space 12} {c |}
{space 4}1.Female {c |}{col 14}{res}{space 2}-1.387598{col 26}{space 2}  .015213{col 37}{space 1}  -91.21{col 46}{space 3}0.000{col 54}{space 4}-1.417454{col 67}{space 3}-1.357741
{txt}std_same_p~x {c |}{col 14}{res}{space 2}-.0221443{col 26}{space 2} .0076728{col 37}{space 1}   -2.89{col 46}{space 3}0.004{col 54}{space 4}-.0372026{col 67}{space 3} -.007086
{txt}{space 12} {c |}
{space 5}AgeGrp4 {c |}
{space 10}2  {c |}{col 14}{res}{space 2} .7791176{col 26}{space 2} .0216724{col 37}{space 1}   35.95{col 46}{space 3}0.000{col 54}{space 4} .7365842{col 67}{space 3} .8216511
{txt}{space 10}3  {c |}{col 14}{res}{space 2} .8912266{col 26}{space 2} .0254944{col 37}{space 1}   34.96{col 46}{space 3}0.000{col 54}{space 4} .8411922{col 67}{space 3}  .941261
{txt}{space 10}4  {c |}{col 14}{res}{space 2} .6804063{col 26}{space 2} .0276824{col 37}{space 1}   24.58{col 46}{space 3}0.000{col 54}{space 4} .6260778{col 67}{space 3} .7347348
{txt}{space 12} {c |}
std_same_p~4 {c |}{col 14}{res}{space 2}-.0389425{col 26}{space 2} .0058551{col 37}{space 1}   -6.65{col 46}{space 3}0.000{col 54}{space 4}-.0504335{col 67}{space 3}-.0274516
{txt}{space 12} {c |}
{space 7}Race5 {c |}
{space 10}2  {c |}{col 14}{res}{space 2}-1.026706{col 26}{space 2} .0309922{col 37}{space 1}  -33.13{col 46}{space 3}0.000{col 54}{space 4} -1.08753{col 67}{space 3}-.9658815
{txt}{space 10}3  {c |}{col 14}{res}{space 2} .0160809{col 26}{space 2} .0791083{col 37}{space 1}    0.20{col 46}{space 3}0.839{col 54}{space 4} -.139174{col 67}{space 3} .1713357
{txt}{space 10}4  {c |}{col 14}{res}{space 2} -.464137{col 26}{space 2} .0535693{col 37}{space 1}   -8.66{col 46}{space 3}0.000{col 54}{space 4}  -.56927{col 67}{space 3} -.359004
{txt}{space 10}5  {c |}{col 14}{res}{space 2}-.7203426{col 26}{space 2} .0417164{col 37}{space 1}  -17.27{col 46}{space 3}0.000{col 54}{space 4}-.8022135{col 67}{space 3}-.6384717
{txt}{space 12} {c |}
std_same_~e5 {c |}{col 14}{res}{space 2} .0128113{col 26}{space 2} .0101364{col 37}{space 1}    1.26{col 46}{space 3}0.207{col 54}{space 4} -.007082{col 67}{space 3} .0327045
{txt}{space 3}1.BornUSA {c |}{col 14}{res}{space 2} .5569525{col 26}{space 2} .0391668{col 37}{space 1}   14.22{col 46}{space 3}0.000{col 54}{space 4} .4800852{col 67}{space 3} .6338197
{txt}std_same_p~A {c |}{col 14}{res}{space 2} .0222712{col 26}{space 2} .0076223{col 37}{space 1}    2.92{col 46}{space 3}0.004{col 54}{space 4} .0073119{col 67}{space 3} .0372305
{txt}{space 12} {c |}
{space 4}MarStat5 {c |}
{space 10}2  {c |}{col 14}{res}{space 2} .8209952{col 26}{space 2} .0214856{col 37}{space 1}   38.21{col 46}{space 3}0.000{col 54}{space 4}  .778828{col 67}{space 3} .8631623
{txt}{space 10}3  {c |}{col 14}{res}{space 2} 1.129914{col 26}{space 2} .0144015{col 37}{space 1}   78.46{col 46}{space 3}0.000{col 54}{space 4}  1.10165{col 67}{space 3} 1.158178
{txt}{space 10}4  {c |}{col 14}{res}{space 2} -.073785{col 26}{space 2} .0905805{col 37}{space 1}   -0.81{col 46}{space 3}0.416{col 54}{space 4}-.2515555{col 67}{space 3} .1039855
{txt}{space 10}5  {c |}{col 14}{res}{space 2} .8250642{col 26}{space 2} .0178972{col 37}{space 1}   46.10{col 46}{space 3}0.000{col 54}{space 4} .7899398{col 67}{space 3} .8601887
{txt}{space 12} {c |}
std_same_~t5 {c |}{col 14}{res}{space 2}-.0316433{col 26}{space 2} .0061531{col 37}{space 1}   -5.14{col 46}{space 3}0.000{col 54}{space 4}-.0437193{col 67}{space 3}-.0195672
{txt}{space 1}Rat_Poverty {c |}{col 14}{res}{space 2} .0036681{col 26}{space 2} .0019263{col 37}{space 1}    1.90{col 46}{space 3}0.057{col 54}{space 4}-.0001124{col 67}{space 3} .0074486
{txt}{space 1}Rat_Mig_Cum {c |}{col 14}{res}{space 2}-1.611869{col 26}{space 2} .4990354{col 37}{space 1}   -3.23{col 46}{space 3}0.001{col 54}{space 4}-2.591256{col 67}{space 3} -.632482
{txt}{space 5}Pop_Den {c |}{col 14}{res}{space 2}-5.62e-06{col 26}{space 2} 4.88e-06{col 37}{space 1}   -1.15{col 46}{space 3}0.250{col 54}{space 4}-.0000152{col 67}{space 3} 3.96e-06
{txt}{space 1}Rat_GC_ProE {c |}{col 14}{res}{space 2} .0000129{col 26}{space 2} 9.14e-06{col 37}{space 1}    1.41{col 46}{space 3}0.160{col 54}{space 4}-5.07e-06{col 67}{space 3} .0000308
{txt}{space 1}Rat_GC_ProM {c |}{col 14}{res}{space 2}-3.73e-07{col 26}{space 2} .0000115{col 37}{space 1}   -0.03{col 46}{space 3}0.974{col 54}{space 4}-.0000228{col 67}{space 3} .0000221
{txt}{space 1}Rat_GC_ProB {c |}{col 14}{res}{space 2}-.0000335{col 26}{space 2}  .000037{col 37}{space 1}   -0.90{col 46}{space 3}0.366{col 54}{space 4}-.0001062{col 67}{space 3} .0000392
{txt}{space 2}Rat_GC_Jew {c |}{col 14}{res}{space 2} .0001161{col 26}{space 2} .0000674{col 37}{space 1}    1.72{col 46}{space 3}0.085{col 54}{space 4}-.0000162{col 67}{space 3} .0002483
{txt}{space 2}Rat_GC_Oth {c |}{col 14}{res}{space 2}-.0001096{col 26}{space 2} .0000231{col 37}{space 1}   -4.74{col 46}{space 3}0.000{col 54}{space 4}-.0001549{col 67}{space 3}-.0000642
{txt}{space 7}_cons {c |}{col 14}{res}{space 2}-9.249288{col 26}{space 2} .1118734{col 37}{space 1}  -82.68{col 46}{space 3}0.000{col 54}{space 4}-9.468847{col 67}{space 3} -9.02973
{txt}{hline 13}{c BT}{hline 11}{hline 11}{hline 9}{hline 8}{hline 13}{hline 12}
{res}{com}.         estimates store m1 
. 
.         if (`model_version' <= 4){c -(}
.                 estimates restore m1
{txt}(results {stata estimates replay m1:m1} are active now)
{com}.                 mchange_mi Female "binary" "mi"

                           {txt}Female
{hline 61}
      Percentiles      Smallest
 1%    {res}        0              0
{txt} 5%    {res}        0              0
{txt}10%    {res}        0              0       {txt}Obs         {res}135,582,733
{txt}25%    {res}        0              0       {txt}Sum of Wgt. {res}  135582733

{txt}50%    {res}        1                      {txt}Mean          {res} .5158124
                        {txt}Largest       Std. Dev.     {res} .4997499
{txt}75%    {res}        1              1
{txt}90%    {res}        1              1       {txt}Variance      {res}   .24975
{txt}95%    {res}        1              1       {txt}Skewness      {res}-.0632814
{txt}99%    {res}        1              1       {txt}Kurtosis      {res} 1.004005

{txt}Multiple-imputation estimates{col 47}Imputations{col 65}= {res}          10
{txt}Predictive margins{col 47}Number of obs{col 65}= {res}  11,764,378

{txt}{col 1}Number of strata{col 19}= {res}        1{txt}{col 47}Population size{col 65}={res}  418,786,443
{txt}{col 1}Number of PSUs{col 19}= {res}      918{txt}{col 47}Subpop. no. obs{col 65}={res}   11,079,265
{txt}{col 47}Subpop. size{col 65}={res}            .
{txt}{col 47}Average RVI{col 65}= {res}      0.0000
{txt}{col 47}Largest FMI{col 65}= {res}      0.0001
{txt}{col 47}Complete DF{col 65}= {res}         917
{txt}DF adjustment:{ralign 15: {res:Small sample}}{col 47}DF:     min{col 65}= {res}      914.94
{txt}{col 47}        avg{col 65}= {res}      914.96
{txt}Within VCE type: {ralign 12:{res:Delta-method}}{col 47}        max{col 65}= {res}      914.98

{txt}Expression{col 14}: {res}Pr(Suic), predict(pr)

{txt}{p2colset 1 14 16 2}{...}
{p2col:1._at}:{space 1}{res:{txt:Female}{space 10}{txt:=} {space 10}0}{p_end}
{p2colreset}{...}

{txt}{p2colset 1 14 16 2}{...}
{p2col:2._at}:{space 1}{res:{txt:Female}{space 10}{txt:=} {space 10}1}{p_end}
{p2colreset}{...}

{res}{txt}{hline 13}{c TT}{hline 11}{hline 11}{hline 9}{hline 8}{hline 13}{hline 12}
{col 14}{c |}     Margin{col 26}   Std. Err.{col 38}      t{col 46}   P>|t|{col 54}     [95% Con{col 67}f. Interval]
{hline 13}{c +}{hline 11}{hline 11}{hline 9}{hline 8}{hline 13}{hline 12}
{space 9}_at {c |}
{space 10}1  {c |}{col 14}{res}{space 2} .0002451{col 26}{space 2} 2.11e-06{col 37}{space 1}  116.40{col 46}{space 3}0.000{col 54}{space 4}  .000241{col 67}{space 3} .0002492
{txt}{space 10}2  {c |}{col 14}{res}{space 2} .0000612{col 26}{space 2} 7.61e-07{col 37}{space 1}   80.45{col 46}{space 3}0.000{col 54}{space 4} .0000597{col 67}{space 3} .0000627
{txt}{hline 13}{c BT}{hline 11}{hline 11}{hline 9}{hline 8}{hline 13}{hline 12}
{res}
{txt}{space 0}{space 0}{ralign 12:}{space 1}{c |}{space 1}{ralign 8:lincom}{space 1}{space 1}{ralign 8:se}{space 1}{space 1}{ralign 8:zvalue}{space 1}{space 1}{ralign 8:pvalue}{space 1}{space 1}{ralign 8:ll}{space 1}{space 1}{ralign 8:ul}{space 1}
{space 0}{hline 13}{c   +}{hline 10}{hline 10}{hline 10}{hline 10}{hline 10}{hline 10}
{space 0}{space 0}{ralign 12:1}{space 1}{c |}{space 1}{ralign 8:{res:{sf:-0.0001839}}}{space 1}{space 1}{ralign 8:{res:{sf:0.0000022}}}{space 1}{space 1}{ralign 8:{res:{sf:-8.2e+01}}}{space 1}{space 1}{ralign 8:{res:{sf: 0.0e+00}}}{space 1}{space 1}{ralign 8:{res:{sf:-0.0001883}}}{space 1}{space 1}{ralign 8:{res:{sf:-0.0001795}}}{space 1}
{com}.                 estimates restore m1
{txt}(results {stata estimates replay m1:m1} are active now)
{com}.                 mchange_mi AgeGrp4 "categorical" "mi"
{res}
{txt}Multiple-imputation estimates{col 47}Imputations{col 65}= {res}          10
{txt}Pairwise comparisons of predictive margins{col 47}Number of obs{col 65}= {res}  11,079,265

{txt}{col 1}Number of strata{col 19}= {res}        1{txt}{col 47}Population size{col 65}={res}  418,786,443
{txt}{col 1}Number of PSUs{col 19}= {res}      918{txt}{col 47}Subpop. no. obs{col 65}={res}   11,079,265
{txt}{col 47}Subpop. size{col 65}={res}            .
{txt}{col 47}Average RVI{col 65}= {res}      0.0005
{txt}{col 47}Largest FMI{col 65}= {res}      0.0011
{txt}{col 47}Complete DF{col 65}= {res}         917
{txt}DF adjustment:{ralign 15: {res:Small sample}}{col 47}DF:     min{col 65}= {res}      913.87
{txt}{col 47}        avg{col 65}= {res}      914.34
{txt}Within VCE type: {ralign 12:{res:Delta-method}}{col 47}        max{col 65}= {res}      914.95

{txt}Expression{col 14}: {res}Pr(Suic), predict(pr)

{txt}{hline 13}{c TT}{hline 11}{hline 11}{hline 14}{hline 12}
{col 14}{c |}   Contrast{col 26}   Std. Err.{col 38}     [95% Con{col 51}f. Interval]
{hline 13}{c +}{hline 11}{hline 11}{hline 14}{hline 12}
{space 5}AgeGrp4 {c |}
{space 5}2 vs 1  {c |}{col 14}{res}{space 2}   .00009{col 26}{space 2} 2.13e-06{col 37}{space 5} .0000859{col 51}{space 3} .0000942
{txt}{space 5}3 vs 1  {c |}{col 14}{res}{space 2} .0001098{col 26}{space 2} 2.77e-06{col 37}{space 5} .0001043{col 51}{space 3} .0001152
{txt}{space 5}4 vs 1  {c |}{col 14}{res}{space 2} .0000744{col 26}{space 2} 2.95e-06{col 37}{space 5} .0000686{col 51}{space 3} .0000802
{txt}{space 5}3 vs 2  {c |}{col 14}{res}{space 2} .0000197{col 26}{space 2} 2.35e-06{col 37}{space 5} .0000151{col 51}{space 3} .0000244
{txt}{space 5}4 vs 2  {c |}{col 14}{res}{space 2}-.0000156{col 26}{space 2} 2.92e-06{col 37}{space 5}-.0000214{col 51}{space 3}-9.90e-06
{txt}{space 5}4 vs 3  {c |}{col 14}{res}{space 2}-.0000354{col 26}{space 2} 2.74e-06{col 37}{space 5}-.0000407{col 51}{space 3}  -.00003
{txt}{hline 13}{c BT}{hline 11}{hline 11}{hline 14}{hline 12}
{res}{com}.                 estimates restore m1
{txt}(results {stata estimates replay m1:m1} are active now)
{com}.                 mchange_mi Race5 "categorical" "mi"
{res}
{txt}Multiple-imputation estimates{col 47}Imputations{col 65}= {res}          10
{txt}Pairwise comparisons of predictive margins{col 47}Number of obs{col 65}= {res}  11,079,265

{txt}{col 1}Number of strata{col 19}= {res}        1{txt}{col 47}Population size{col 65}={res}  418,786,443
{txt}{col 1}Number of PSUs{col 19}= {res}      918{txt}{col 47}Subpop. no. obs{col 65}={res}   11,079,265
{txt}{col 47}Subpop. size{col 65}={res}            .
{txt}{col 47}Average RVI{col 65}= {res}      0.0000
{txt}{col 47}Largest FMI{col 65}= {res}      0.0000
{txt}{col 47}Complete DF{col 65}= {res}         917
{txt}DF adjustment:{ralign 15: {res:Small sample}}{col 47}DF:     min{col 65}= {res}      914.98
{txt}{col 47}        avg{col 65}= {res}      915.00
{txt}Within VCE type: {ralign 12:{res:Delta-method}}{col 47}        max{col 65}= {res}      915.01

{txt}Expression{col 14}: {res}Pr(Suic), predict(pr)

{txt}{hline 13}{c TT}{hline 11}{hline 11}{hline 14}{hline 12}
{col 14}{c |}   Contrast{col 26}   Std. Err.{col 38}     [95% Con{col 51}f. Interval]
{hline 13}{c +}{hline 11}{hline 11}{hline 14}{hline 12}
{space 7}Race5 {c |}
{space 5}2 vs 1  {c |}{col 14}{res}{space 2}-.0001117{col 26}{space 2} 2.66e-06{col 37}{space 5}-.0001169{col 51}{space 3}-.0001065
{txt}{space 5}3 vs 1  {c |}{col 14}{res}{space 2} 2.82e-06{col 26}{space 2}  .000014{col 37}{space 5}-.0000246{col 51}{space 3} .0000303
{txt}{space 5}4 vs 1  {c |}{col 14}{res}{space 2}-.0000646{col 26}{space 2} 6.16e-06{col 37}{space 5}-.0000767{col 51}{space 3}-.0000525
{txt}{space 5}5 vs 1  {c |}{col 14}{res}{space 2}-.0000893{col 26}{space 2} 3.94e-06{col 37}{space 5}-.0000971{col 51}{space 3}-.0000816
{txt}{space 5}3 vs 2  {c |}{col 14}{res}{space 2} .0001145{col 26}{space 2} .0000146{col 37}{space 5} .0000859{col 51}{space 3} .0001431
{txt}{space 5}4 vs 2  {c |}{col 14}{res}{space 2} .0000471{col 26}{space 2} 5.17e-06{col 37}{space 5} .0000369{col 51}{space 3} .0000572
{txt}{space 5}5 vs 2  {c |}{col 14}{res}{space 2} .0000223{col 26}{space 2} 3.37e-06{col 37}{space 5} .0000157{col 51}{space 3}  .000029
{txt}{space 5}4 vs 3  {c |}{col 14}{res}{space 2}-.0000674{col 26}{space 2} .0000163{col 37}{space 5}-.0000993{col 51}{space 3}-.0000355
{txt}{space 5}5 vs 3  {c |}{col 14}{res}{space 2}-.0000922{col 26}{space 2} .0000153{col 37}{space 5}-.0001222{col 51}{space 3}-.0000621
{txt}{space 5}5 vs 4  {c |}{col 14}{res}{space 2}-.0000247{col 26}{space 2} 5.41e-06{col 37}{space 5}-.0000353{col 51}{space 3}-.0000141
{txt}{hline 13}{c BT}{hline 11}{hline 11}{hline 14}{hline 12}
{res}{com}.                 estimates restore m1
{txt}(results {stata estimates replay m1:m1} are active now)
{com}.                 mchange_mi BornUSA "binary" "mi"

                           {txt}BornUSA
{hline 61}
      Percentiles      Smallest
 1%    {res}        0              0
{txt} 5%    {res}        0              0
{txt}10%    {res}        1              0       {txt}Obs         {res}135,582,733
{txt}25%    {res}        1              0       {txt}Sum of Wgt. {res}  135582733

{txt}50%    {res}        1                      {txt}Mean          {res} .9265123
                        {txt}Largest       Std. Dev.     {res} .2609353
{txt}75%    {res}        1              1
{txt}90%    {res}        1              1       {txt}Variance      {res} .0680872
{txt}95%    {res}        1              1       {txt}Skewness      {res}-3.269105
{txt}99%    {res}        1              1       {txt}Kurtosis      {res} 11.68704

{txt}Multiple-imputation estimates{col 47}Imputations{col 65}= {res}          10
{txt}Predictive margins{col 47}Number of obs{col 65}= {res}  11,764,378

{txt}{col 1}Number of strata{col 19}= {res}        1{txt}{col 47}Population size{col 65}={res}  418,786,443
{txt}{col 1}Number of PSUs{col 19}= {res}      918{txt}{col 47}Subpop. no. obs{col 65}={res}   11,079,265
{txt}{col 47}Subpop. size{col 65}={res}            .
{txt}{col 47}Average RVI{col 65}= {res}      0.0000
{txt}{col 47}Largest FMI{col 65}= {res}      0.0000
{txt}{col 47}Complete DF{col 65}= {res}         917
{txt}DF adjustment:{ralign 15: {res:Small sample}}{col 47}DF:     min{col 65}= {res}      915.00
{txt}{col 47}        avg{col 65}= {res}      915.00
{txt}Within VCE type: {ralign 12:{res:Delta-method}}{col 47}        max{col 65}= {res}      915.01

{txt}Expression{col 14}: {res}Pr(Suic), predict(pr)

{txt}{p2colset 1 14 16 2}{...}
{p2col:1._at}:{space 1}{res:{txt:BornUSA}{space 9}{txt:=} {space 10}0}{p_end}
{p2colreset}{...}

{txt}{p2colset 1 14 16 2}{...}
{p2col:2._at}:{space 1}{res:{txt:BornUSA}{space 9}{txt:=} {space 10}1}{p_end}
{p2colreset}{...}

{res}{txt}{hline 13}{c TT}{hline 11}{hline 11}{hline 9}{hline 8}{hline 13}{hline 12}
{col 14}{c |}     Margin{col 26}   Std. Err.{col 38}      t{col 46}   P>|t|{col 54}     [95% Con{col 67}f. Interval]
{hline 13}{c +}{hline 11}{hline 11}{hline 9}{hline 8}{hline 13}{hline 12}
{space 9}_at {c |}
{space 10}1  {c |}{col 14}{res}{space 2} .0000886{col 26}{space 2} 3.17e-06{col 37}{space 1}   27.91{col 46}{space 3}0.000{col 54}{space 4} .0000824{col 67}{space 3} .0000948
{txt}{space 10}2  {c |}{col 14}{res}{space 2} .0001546{col 26}{space 2} 1.24e-06{col 37}{space 1}  124.40{col 46}{space 3}0.000{col 54}{space 4} .0001521{col 67}{space 3}  .000157
{txt}{hline 13}{c BT}{hline 11}{hline 11}{hline 9}{hline 8}{hline 13}{hline 12}
{res}
{txt}{space 0}{space 0}{ralign 12:}{space 1}{c |}{space 1}{ralign 8:lincom}{space 1}{space 1}{ralign 8:se}{space 1}{space 1}{ralign 8:zvalue}{space 1}{space 1}{ralign 8:pvalue}{space 1}{space 1}{ralign 8:ll}{space 1}{space 1}{ralign 8:ul}{space 1}
{space 0}{hline 13}{c   +}{hline 10}{hline 10}{hline 10}{hline 10}{hline 10}{hline 10}
{space 0}{space 0}{ralign 12:1}{space 1}{c |}{space 1}{ralign 8:{res:{sf:0.0000660}}}{space 1}{space 1}{ralign 8:{res:{sf:0.0000038}}}{space 1}{space 1}{ralign 8:{res:{sf: 1.8e+01}}}{space 1}{space 1}{ralign 8:{res:{sf: 0.0e+00}}}{space 1}{space 1}{ralign 8:{res:{sf:0.0000586}}}{space 1}{space 1}{ralign 8:{res:{sf:0.0000734}}}{space 1}
{com}.                 estimates restore m1
{txt}(results {stata estimates replay m1:m1} are active now)
{com}.                 mchange_mi MarStat5 "categorical" "mi"
{res}
{txt}Multiple-imputation estimates{col 47}Imputations{col 65}= {res}          10
{txt}Pairwise comparisons of predictive margins{col 47}Number of obs{col 65}= {res}  11,079,265

{txt}{col 1}Number of strata{col 19}= {res}        1{txt}{col 47}Population size{col 65}={res}  418,786,443
{txt}{col 1}Number of PSUs{col 19}= {res}      918{txt}{col 47}Subpop. no. obs{col 65}={res}   11,079,265
{txt}{col 47}Subpop. size{col 65}={res}            .
{txt}{col 47}Average RVI{col 65}= {res}      0.0030
{txt}{col 47}Largest FMI{col 65}= {res}      0.0038
{txt}{col 47}Complete DF{col 65}= {res}         917
{txt}DF adjustment:{ralign 15: {res:Small sample}}{col 47}DF:     min{col 65}= {res}      910.20
{txt}{col 47}        avg{col 65}= {res}      911.51
{txt}Within VCE type: {ralign 12:{res:Delta-method}}{col 47}        max{col 65}= {res}      913.43

{txt}Expression{col 14}: {res}Pr(Suic), predict(pr)

{txt}{hline 13}{c TT}{hline 11}{hline 11}{hline 14}{hline 12}
{col 14}{c |}   Contrast{col 26}   Std. Err.{col 38}     [95% Con{col 51}f. Interval]
{hline 13}{c +}{hline 11}{hline 11}{hline 14}{hline 12}
{space 4}MarStat5 {c |}
{space 5}2 vs 1  {c |}{col 14}{res}{space 2} .0001193{col 26}{space 2} 4.36e-06{col 37}{space 5} .0001107{col 51}{space 3} .0001279
{txt}{space 5}3 vs 1  {c |}{col 14}{res}{space 2} .0001964{col 26}{space 2} 3.30e-06{col 37}{space 5} .0001899{col 51}{space 3} .0002029
{txt}{space 5}4 vs 1  {c |}{col 14}{res}{space 2}-6.67e-06{col 26}{space 2} 7.91e-06{col 37}{space 5}-.0000222{col 51}{space 3} 8.85e-06
{txt}{space 5}5 vs 1  {c |}{col 14}{res}{space 2} .0001202{col 26}{space 2} 2.88e-06{col 37}{space 5} .0001145{col 51}{space 3} .0001258
{txt}{space 5}3 vs 2  {c |}{col 14}{res}{space 2} .0000771{col 26}{space 2} 5.06e-06{col 37}{space 5} .0000672{col 51}{space 3}  .000087
{txt}{space 5}4 vs 2  {c |}{col 14}{res}{space 2} -.000126{col 26}{space 2} 9.07e-06{col 37}{space 5}-.0001438{col 51}{space 3}-.0001082
{txt}{space 5}5 vs 2  {c |}{col 14}{res}{space 2} 8.68e-07{col 26}{space 2} 5.38e-06{col 37}{space 5}-9.69e-06{col 51}{space 3} .0000114
{txt}{space 5}4 vs 3  {c |}{col 14}{res}{space 2}-.0002031{col 26}{space 2} 8.12e-06{col 37}{space 5} -.000219{col 51}{space 3}-.0001871
{txt}{space 5}5 vs 3  {c |}{col 14}{res}{space 2}-.0000762{col 26}{space 2} 4.34e-06{col 37}{space 5}-.0000847{col 51}{space 3}-.0000677
{txt}{space 5}5 vs 4  {c |}{col 14}{res}{space 2} .0001268{col 26}{space 2} 7.77e-06{col 37}{space 5} .0001116{col 51}{space 3} .0001421
{txt}{hline 13}{c BT}{hline 11}{hline 11}{hline 14}{hline 12}
{res}{com}.                 
.                 if (`model_version' == 1 | `model_version' == 2) {c -(}
.                         estimates restore m1
.                         mchange_mi RAT_Female "continuous" "mi"
.                         estimates restore m1
.                         mchange_mi RAT_AgeGrp4_2 "continuous" "mi"
.                         estimates restore m1
.                         mchange_mi RAT_AgeGrp4_3 "continuous" "mi"
.                         estimates restore m1
.                         mchange_mi RAT_AgeGrp4_4 "continuous" "mi"
.                         estimates restore m1
.                         mchange_mi RAT_Race5_2 "continuous" "mi"
.                         estimates restore m1
.                         mchange_mi RAT_Race5_3 "continuous" "mi"
.                         estimates restore m1
.                         mchange_mi RAT_Race5_4 "continuous" "mi"
.                         estimates restore m1
.                         mchange_mi RAT_Race5_5 "continuous" "mi"
.                         estimates restore m1
.                         mchange_mi RAT_BornUSA "continuous" "mi"
.                         estimates restore m1
.                         mchange_mi RAT_MarStat5_2 "continuous" "mi"
.                         estimates restore m1
.                         mchange_mi RAT_MarStat5_3 "continuous" "mi"
.                         estimates restore m1
.                         mchange_mi RAT_MarStat5_4 "continuous" "mi"
.                         estimates restore m1
.                         mchange_mi RAT_MarStat5_5 "continuous" "mi"
.                 {c )-}
.                 if (`model_version' == 3 | `model_version' == 4) {c -(}
.                         estimates restore m1
{txt}(results {stata estimates replay m1:m1} are active now)
{com}.                         mchange_mi std_same_prop_Sex "continuous" "mi"

                      {txt}std_same_prop_Sex
{hline 61}
      Percentiles      Smallest
 1%    {res}-2.509072      -4.036727
{txt} 5%    {res}-1.338853      -4.036727
{txt}10%    {res}-.9852811      -4.036727       {txt}Obs         {res}135,582,733
{txt}25%    {res}-.5357422      -4.036727       {txt}Sum of Wgt. {res}  135582733

{txt}50%    {res} .0312444                      {txt}Mean          {res}  .027466
                        {txt}Largest       Std. Dev.     {res} .9119966
{txt}75%    {res} .5763232       4.036727
{txt}90%    {res} 1.032018       4.036727       {txt}Variance      {res} .8317377
{txt}95%    {res} 1.401881       4.036727       {txt}Skewness      {res} .0293361
{txt}99%    {res} 2.682605       4.036727       {txt}Kurtosis      {res} 4.966479

{txt}Multiple-imputation estimates{col 47}Imputations{col 65}= {res}          10
{txt}Predictive margins{col 47}Number of obs{col 65}= {res}  11,764,378

{txt}{col 1}Number of strata{col 19}= {res}        1{txt}{col 47}Population size{col 65}={res}  418,786,443
{txt}{col 1}Number of PSUs{col 19}= {res}      918{txt}{col 47}Subpop. no. obs{col 65}={res}   11,079,265
{txt}{col 47}Subpop. size{col 65}={res}            .
{txt}{col 47}Average RVI{col 65}= {res}      0.0000
{txt}{col 47}Largest FMI{col 65}= {res}      0.0000
{txt}{col 47}Complete DF{col 65}= {res}         917
{txt}DF adjustment:{ralign 15: {res:Small sample}}{col 47}DF:     min{col 65}= {res}      915.01
{txt}{col 47}        avg{col 65}= {res}      915.01
{txt}Within VCE type: {ralign 12:{res:Delta-method}}{col 47}        max{col 65}= {res}      915.01

{txt}Expression{col 14}: {res}Pr(Suic), predict(pr)

{txt}{p2colset 1 14 16 2}{...}
{p2col:1._at}:{space 1}{res:{txt:std_same_p~x}{space 4}{txt:=} {space 2}-4.036727}{p_end}
{p2colreset}{...}

{txt}{p2colset 1 14 16 2}{...}
{p2col:2._at}:{space 1}{res:{txt:std_same_p~x}{space 4}{txt:=} {space 3}4.036727}{p_end}
{p2colreset}{...}

{res}{txt}{hline 13}{c TT}{hline 11}{hline 11}{hline 9}{hline 8}{hline 13}{hline 12}
{col 14}{c |}     Margin{col 26}   Std. Err.{col 38}      t{col 46}   P>|t|{col 54}     [95% Con{col 67}f. Interval]
{hline 13}{c +}{hline 11}{hline 11}{hline 9}{hline 8}{hline 13}{hline 12}
{space 9}_at {c |}
{space 10}1  {c |}{col 14}{res}{space 2} .0001619{col 26}{space 2} 5.05e-06{col 37}{space 1}   32.06{col 46}{space 3}0.000{col 54}{space 4}  .000152{col 67}{space 3} .0001718
{txt}{space 10}2  {c |}{col 14}{res}{space 2} .0001354{col 26}{space 2} 4.38e-06{col 37}{space 1}   30.91{col 46}{space 3}0.000{col 54}{space 4} .0001268{col 67}{space 3}  .000144
{txt}{hline 13}{c BT}{hline 11}{hline 11}{hline 9}{hline 8}{hline 13}{hline 12}
{res}
{txt}{space 0}{space 0}{ralign 12:}{space 1}{c |}{space 1}{ralign 8:lincom}{space 1}{space 1}{ralign 8:se}{space 1}{space 1}{ralign 8:zvalue}{space 1}{space 1}{ralign 8:pvalue}{space 1}{space 1}{ralign 8:ll}{space 1}{space 1}{ralign 8:ul}{space 1}
{space 0}{hline 13}{c   +}{hline 10}{hline 10}{hline 10}{hline 10}{hline 10}{hline 10}
{space 0}{space 0}{ralign 12:1}{space 1}{c |}{space 1}{ralign 8:{res:{sf:-0.0000265}}}{space 1}{space 1}{ralign 8:{res:{sf:0.0000092}}}{space 1}{space 1}{ralign 8:{res:{sf:-2.9e+00}}}{space 1}{space 1}{ralign 8:{res:{sf:0.0040344}}}{space 1}{space 1}{ralign 8:{res:{sf:-0.0000445}}}{space 1}{space 1}{ralign 8:{res:{sf:-0.0000085}}}{space 1}
{com}.                         estimates restore m1
{txt}(results {stata estimates replay m1:m1} are active now)
{com}.                         mchange_mi std_same_prop_AgeGrp4 "continuous" "mi"

                    {txt}std_same_prop_AgeGrp4
{hline 61}
      Percentiles      Smallest
 1%    {res}-2.469153      -5.047376
{txt} 5%    {res}-1.445343      -5.047376
{txt}10%    {res}-1.068901      -5.047376       {txt}Obs         {res}135,582,733
{txt}25%    {res} -.504598      -5.047376       {txt}Sum of Wgt. {res}  135582733

{txt}50%    {res} .0808892                      {txt}Mean          {res}  .127631
                        {txt}Largest       Std. Dev.     {res} 1.053953
{txt}75%    {res} .6954771       6.771116
{txt}90%    {res} 1.347035       6.771116       {txt}Variance      {res} 1.110817
{txt}95%    {res} 1.891168       6.771116       {txt}Skewness      {res} .3782019
{txt}99%    {res} 2.953919       6.771116       {txt}Kurtosis      {res}  5.03041

{txt}Multiple-imputation estimates{col 47}Imputations{col 65}= {res}          10
{txt}Predictive margins{col 47}Number of obs{col 65}= {res}  11,764,378

{txt}{col 1}Number of strata{col 19}= {res}        1{txt}{col 47}Population size{col 65}={res}  418,786,443
{txt}{col 1}Number of PSUs{col 19}= {res}      918{txt}{col 47}Subpop. no. obs{col 65}={res}   11,079,265
{txt}{col 47}Subpop. size{col 65}={res}            .
{txt}{col 47}Average RVI{col 65}= {res}      0.0001
{txt}{col 47}Largest FMI{col 65}= {res}      0.0001
{txt}{col 47}Complete DF{col 65}= {res}         917
{txt}DF adjustment:{ralign 15: {res:Small sample}}{col 47}DF:     min{col 65}= {res}      914.88
{txt}{col 47}        avg{col 65}= {res}      914.88
{txt}Within VCE type: {ralign 12:{res:Delta-method}}{col 47}        max{col 65}= {res}      914.89

{txt}Expression{col 14}: {res}Pr(Suic), predict(pr)

{txt}{p2colset 1 14 16 2}{...}
{p2col:1._at}:{space 1}{res:{txt:std_same_p~4}{space 4}{txt:=} {space 2}-5.047376}{p_end}
{p2colreset}{...}

{txt}{p2colset 1 14 16 2}{...}
{p2col:2._at}:{space 1}{res:{txt:std_same_p~4}{space 4}{txt:=} {space 3}6.771116}{p_end}
{p2colreset}{...}

{res}{txt}{hline 13}{c TT}{hline 11}{hline 11}{hline 9}{hline 8}{hline 13}{hline 12}
{col 14}{c |}     Margin{col 26}   Std. Err.{col 38}      t{col 46}   P>|t|{col 54}     [95% Con{col 67}f. Interval]
{hline 13}{c +}{hline 11}{hline 11}{hline 9}{hline 8}{hline 13}{hline 12}
{space 9}_at {c |}
{space 10}1  {c |}{col 14}{res}{space 2} .0001819{col 26}{space 2} 5.50e-06{col 37}{space 1}   33.09{col 46}{space 3}0.000{col 54}{space 4} .0001711{col 67}{space 3} .0001927
{txt}{space 10}2  {c |}{col 14}{res}{space 2} .0001148{col 26}{space 2} 4.64e-06{col 37}{space 1}   24.74{col 46}{space 3}0.000{col 54}{space 4} .0001057{col 67}{space 3} .0001239
{txt}{hline 13}{c BT}{hline 11}{hline 11}{hline 9}{hline 8}{hline 13}{hline 12}
{res}
{txt}{space 0}{space 0}{ralign 12:}{space 1}{c |}{space 1}{ralign 8:lincom}{space 1}{space 1}{ralign 8:se}{space 1}{space 1}{ralign 8:zvalue}{space 1}{space 1}{ralign 8:pvalue}{space 1}{space 1}{ralign 8:ll}{space 1}{space 1}{ralign 8:ul}{space 1}
{space 0}{hline 13}{c   +}{hline 10}{hline 10}{hline 10}{hline 10}{hline 10}{hline 10}
{space 0}{space 0}{ralign 12:1}{space 1}{c |}{space 1}{ralign 8:{res:{sf:-0.0000671}}}{space 1}{space 1}{ralign 8:{res:{sf:0.0000099}}}{space 1}{space 1}{ralign 8:{res:{sf:-6.8e+00}}}{space 1}{space 1}{ralign 8:{res:{sf:0.0000000}}}{space 1}{space 1}{ralign 8:{res:{sf:-0.0000866}}}{space 1}{space 1}{ralign 8:{res:{sf:-0.0000476}}}{space 1}
{com}.                         estimates restore m1
{txt}(results {stata estimates replay m1:m1} are active now)
{com}.                         mchange_mi std_same_prop_Race5 "continuous" "mi"

                     {txt}std_same_prop_Race5
{hline 61}
      Percentiles      Smallest
 1%    {res}-2.270736      -3.614547
{txt} 5%    {res}-1.387242      -3.614547
{txt}10%    {res}-.9511417      -3.614547       {txt}Obs         {res}135,582,733
{txt}25%    {res}-.2595515      -3.614547       {txt}Sum of Wgt. {res}  135582733

{txt}50%    {res} .4283832                      {txt}Mean          {res} .4242603
                        {txt}Largest       Std. Dev.     {res} 1.233905
{txt}75%    {res} 1.013094       8.750521
{txt}90%    {res} 1.238259       8.750521       {txt}Variance      {res} 1.522523
{txt}95%    {res} 2.074838       8.750521       {txt}Skewness      {res} 1.661957
{txt}99%    {res} 6.478749       8.750521       {txt}Kurtosis      {res} 11.30467

{txt}Multiple-imputation estimates{col 47}Imputations{col 65}= {res}          10
{txt}Predictive margins{col 47}Number of obs{col 65}= {res}  11,764,378

{txt}{col 1}Number of strata{col 19}= {res}        1{txt}{col 47}Population size{col 65}={res}  418,786,443
{txt}{col 1}Number of PSUs{col 19}= {res}      918{txt}{col 47}Subpop. no. obs{col 65}={res}   11,079,265
{txt}{col 47}Subpop. size{col 65}={res}            .
{txt}{col 47}Average RVI{col 65}= {res}      0.0000
{txt}{col 47}Largest FMI{col 65}= {res}      0.0000
{txt}{col 47}Complete DF{col 65}= {res}         917
{txt}DF adjustment:{ralign 15: {res:Small sample}}{col 47}DF:     min{col 65}= {res}      915.00
{txt}{col 47}        avg{col 65}= {res}      915.00
{txt}Within VCE type: {ralign 12:{res:Delta-method}}{col 47}        max{col 65}= {res}      915.00

{txt}Expression{col 14}: {res}Pr(Suic), predict(pr)

{txt}{p2colset 1 14 16 2}{...}
{p2col:1._at}:{space 1}{res:{txt:std_same_~e5}{space 4}{txt:=} {space 2}-3.614547}{p_end}
{p2colreset}{...}

{txt}{p2colset 1 14 16 2}{...}
{p2col:2._at}:{space 1}{res:{txt:std_same_~e5}{space 4}{txt:=} {space 3}8.750521}{p_end}
{p2colreset}{...}

{res}{txt}{hline 13}{c TT}{hline 11}{hline 11}{hline 9}{hline 8}{hline 13}{hline 12}
{col 14}{c |}     Margin{col 26}   Std. Err.{col 38}      t{col 46}   P>|t|{col 54}     [95% Con{col 67}f. Interval]
{hline 13}{c +}{hline 11}{hline 11}{hline 9}{hline 8}{hline 13}{hline 12}
{space 9}_at {c |}
{space 10}1  {c |}{col 14}{res}{space 2} .0001414{col 26}{space 2} 5.46e-06{col 37}{space 1}   25.89{col 46}{space 3}0.000{col 54}{space 4} .0001307{col 67}{space 3} .0001522
{txt}{space 10}2  {c |}{col 14}{res}{space 2} .0001657{col 26}{space 2} .0000145{col 37}{space 1}   11.41{col 46}{space 3}0.000{col 54}{space 4} .0001372{col 67}{space 3} .0001942
{txt}{hline 13}{c BT}{hline 11}{hline 11}{hline 9}{hline 8}{hline 13}{hline 12}
{res}
{txt}{space 0}{space 0}{ralign 12:}{space 1}{c |}{space 1}{ralign 8:lincom}{space 1}{space 1}{ralign 8:se}{space 1}{space 1}{ralign 8:zvalue}{space 1}{space 1}{ralign 8:pvalue}{space 1}{space 1}{ralign 8:ll}{space 1}{space 1}{ralign 8:ul}{space 1}
{space 0}{hline 13}{c   +}{hline 10}{hline 10}{hline 10}{hline 10}{hline 10}{hline 10}
{space 0}{space 0}{ralign 12:1}{space 1}{c |}{space 1}{ralign 8:{res:{sf:0.0000243}}}{space 1}{space 1}{ralign 8:{res:{sf:0.0000198}}}{space 1}{space 1}{ralign 8:{res:{sf: 1.2e+00}}}{space 1}{space 1}{ralign 8:{res:{sf:0.2215907}}}{space 1}{space 1}{ralign 8:{res:{sf:-0.0000147}}}{space 1}{space 1}{ralign 8:{res:{sf:0.0000632}}}{space 1}
{com}.                         estimates restore m1
{txt}(results {stata estimates replay m1:m1} are active now)
{com}.                         mchange_mi std_same_prop_BornUSA "continuous" "mi"

                    {txt}std_same_prop_BornUSA
{hline 61}
      Percentiles      Smallest
 1%    {res}-5.045306      -7.936538
{txt} 5%    {res}-2.723737      -7.936538
{txt}10%    {res}-1.754185      -7.936538       {txt}Obs         {res}135,582,733
{txt}25%    {res}-.5298624      -7.936538       {txt}Sum of Wgt. {res}  135582733

{txt}50%    {res} .1534955                      {txt}Mean          {res}-.1095675
                        {txt}Largest       Std. Dev.     {res} 1.400314
{txt}75%    {res} .6318931       7.936538
{txt}90%    {res} .8333296       7.936538       {txt}Variance      {res} 1.960879
{txt}95%    {res} .9703126       7.936538       {txt}Skewness      {res}-.7379761
{txt}99%    {res} 4.547946       7.936538       {txt}Kurtosis      {res} 8.840372

{txt}Multiple-imputation estimates{col 47}Imputations{col 65}= {res}          10
{txt}Predictive margins{col 47}Number of obs{col 65}= {res}  11,764,378

{txt}{col 1}Number of strata{col 19}= {res}        1{txt}{col 47}Population size{col 65}={res}  418,786,443
{txt}{col 1}Number of PSUs{col 19}= {res}      918{txt}{col 47}Subpop. no. obs{col 65}={res}   11,079,265
{txt}{col 47}Subpop. size{col 65}={res}            .
{txt}{col 47}Average RVI{col 65}= {res}      0.0000
{txt}{col 47}Largest FMI{col 65}= {res}      0.0000
{txt}{col 47}Complete DF{col 65}= {res}         917
{txt}DF adjustment:{ralign 15: {res:Small sample}}{col 47}DF:     min{col 65}= {res}      915.01
{txt}{col 47}        avg{col 65}= {res}      915.01
{txt}Within VCE type: {ralign 12:{res:Delta-method}}{col 47}        max{col 65}= {res}      915.01

{txt}Expression{col 14}: {res}Pr(Suic), predict(pr)

{txt}{p2colset 1 14 16 2}{...}
{p2col:1._at}:{space 1}{res:{txt:std_same_p~A}{space 4}{txt:=} {space 2}-7.936538}{p_end}
{p2colreset}{...}

{txt}{p2colset 1 14 16 2}{...}
{p2col:2._at}:{space 1}{res:{txt:std_same_p~A}{space 4}{txt:=} {space 3}7.936538}{p_end}
{p2colreset}{...}

{res}{txt}{hline 13}{c TT}{hline 11}{hline 11}{hline 9}{hline 8}{hline 13}{hline 12}
{col 14}{c |}     Margin{col 26}   Std. Err.{col 38}      t{col 46}   P>|t|{col 54}     [95% Con{col 67}f. Interval]
{hline 13}{c +}{hline 11}{hline 11}{hline 9}{hline 8}{hline 13}{hline 12}
{space 9}_at {c |}
{space 10}1  {c |}{col 14}{res}{space 2}  .000126{col 26}{space 2} 7.19e-06{col 37}{space 1}   17.52{col 46}{space 3}0.000{col 54}{space 4} .0001119{col 67}{space 3} .0001401
{txt}{space 10}2  {c |}{col 14}{res}{space 2} .0001794{col 26}{space 2} .0000116{col 37}{space 1}   15.45{col 46}{space 3}0.000{col 54}{space 4} .0001566{col 67}{space 3} .0002022
{txt}{hline 13}{c BT}{hline 11}{hline 11}{hline 9}{hline 8}{hline 13}{hline 12}
{res}
{txt}{space 0}{space 0}{ralign 12:}{space 1}{c |}{space 1}{ralign 8:lincom}{space 1}{space 1}{ralign 8:se}{space 1}{space 1}{ralign 8:zvalue}{space 1}{space 1}{ralign 8:pvalue}{space 1}{space 1}{ralign 8:ll}{space 1}{space 1}{ralign 8:ul}{space 1}
{space 0}{hline 13}{c   +}{hline 10}{hline 10}{hline 10}{hline 10}{hline 10}{hline 10}
{space 0}{space 0}{ralign 12:1}{space 1}{c |}{space 1}{ralign 8:{res:{sf:0.0000534}}}{space 1}{space 1}{ralign 8:{res:{sf:0.0000187}}}{space 1}{space 1}{ralign 8:{res:{sf: 2.9e+00}}}{space 1}{space 1}{ralign 8:{res:{sf:0.0043423}}}{space 1}{space 1}{ralign 8:{res:{sf:0.0000168}}}{space 1}{space 1}{ralign 8:{res:{sf:0.0000901}}}{space 1}
{com}.                         estimates restore m1
{txt}(results {stata estimates replay m1:m1} are active now)
{com}.                         mchange_mi std_same_prop_MarStat5 "continuous" "mi"

                   {txt}std_same_prop_MarStat5
{hline 61}
      Percentiles      Smallest
 1%    {res}-2.292691      -4.824302
{txt} 5%    {res}-1.574807      -4.824302
{txt}10%    {res}-1.168543      -4.824302       {txt}Obs         {res}135,582,272
{txt}25%    {res}-.5012236      -4.824302       {txt}Sum of Wgt. {res}  135582272

{txt}50%    {res} .1833285                      {txt}Mean          {res} .1423099
                        {txt}Largest       Std. Dev.     {res} 1.051846
{txt}75%    {res} .8066256         5.1647
{txt}90%    {res} 1.365528         5.1647       {txt}Variance      {res}  1.10638
{txt}95%    {res} 1.765109         5.1647       {txt}Skewness      {res} .0184338
{txt}99%    {res} 2.683285         5.1647       {txt}Kurtosis      {res}  4.57374

{txt}Multiple-imputation estimates{col 47}Imputations{col 65}= {res}          10
{txt}Predictive margins{col 47}Number of obs{col 65}= {res}  11,764,378

{txt}{col 1}Number of strata{col 19}= {res}        1{txt}{col 47}Population size{col 65}={res}  418,786,443
{txt}{col 1}Number of PSUs{col 19}= {res}      918{txt}{col 47}Subpop. no. obs{col 65}={res}   11,079,265
{txt}{col 47}Subpop. size{col 65}={res}            .
{txt}{col 47}Average RVI{col 65}= {res}      0.0022
{txt}{col 47}Largest FMI{col 65}= {res}      0.0041
{txt}{col 47}Complete DF{col 65}= {res}         917
{txt}DF adjustment:{ralign 15: {res:Small sample}}{col 47}DF:     min{col 65}= {res}      909.66
{txt}{col 47}        avg{col 65}= {res}      909.72
{txt}Within VCE type: {ralign 12:{res:Delta-method}}{col 47}        max{col 65}= {res}      909.78

{txt}Expression{col 14}: {res}Pr(Suic), predict(pr)

{txt}{p2colset 1 14 16 2}{...}
{p2col:1._at}:{space 1}{res:{txt:std_same_~t5}{space 4}{txt:=} {space 2}-4.824302}{p_end}
{p2colreset}{...}

{txt}{p2colset 1 14 16 2}{...}
{p2col:2._at}:{space 1}{res:{txt:std_same_~t5}{space 4}{txt:=} {space 5}5.1647}{p_end}
{p2colreset}{...}

{res}{txt}{hline 13}{c TT}{hline 11}{hline 11}{hline 9}{hline 8}{hline 13}{hline 12}
{col 14}{c |}     Margin{col 26}   Std. Err.{col 38}      t{col 46}   P>|t|{col 54}     [95% Con{col 67}f. Interval]
{hline 13}{c +}{hline 11}{hline 11}{hline 9}{hline 8}{hline 13}{hline 12}
{space 9}_at {c |}
{space 10}1  {c |}{col 14}{res}{space 2} .0001738{col 26}{space 2} 5.48e-06{col 37}{space 1}   31.73{col 46}{space 3}0.000{col 54}{space 4} .0001631{col 67}{space 3} .0001846
{txt}{space 10}2  {c |}{col 14}{res}{space 2} .0001267{col 26}{space 2} 4.00e-06{col 37}{space 1}   31.68{col 46}{space 3}0.000{col 54}{space 4} .0001189{col 67}{space 3} .0001346
{txt}{hline 13}{c BT}{hline 11}{hline 11}{hline 9}{hline 8}{hline 13}{hline 12}
{res}
{txt}{space 0}{space 0}{ralign 12:}{space 1}{c |}{space 1}{ralign 8:lincom}{space 1}{space 1}{ralign 8:se}{space 1}{space 1}{ralign 8:zvalue}{space 1}{space 1}{ralign 8:pvalue}{space 1}{space 1}{ralign 8:ll}{space 1}{space 1}{ralign 8:ul}{space 1}
{space 0}{hline 13}{c   +}{hline 10}{hline 10}{hline 10}{hline 10}{hline 10}{hline 10}
{space 0}{space 0}{ralign 12:1}{space 1}{c |}{space 1}{ralign 8:{res:{sf:-0.0000471}}}{space 1}{space 1}{ralign 8:{res:{sf:0.0000092}}}{space 1}{space 1}{ralign 8:{res:{sf:-5.1e+00}}}{space 1}{space 1}{ralign 8:{res:{sf:0.0000004}}}{space 1}{space 1}{ralign 8:{res:{sf:-0.0000652}}}{space 1}{space 1}{ralign 8:{res:{sf:-0.0000290}}}{space 1}
{com}.                 {c )-}
.         
.                 if (`model_version' == 2 | `model_version' == 4){c -(}
.                         estimates restore m1
.                         mchange_mi UnEmpl "categorical" "mi"
.                         estimates restore m1
.                         mchange_mi PhysProb "categorical" "mi"
.                 {c )-}
.         
.                 if (`model_version' == 2){c -(}
.                         estimates restore m1
.                         mchange_mi RAT_UnEmpl "continuous" "mi"
.                         estimates restore m1
.                         mchange_mi RAT_PhysProb "continuous" "mi"
.                 {c )-}
.                 if (`model_version' == 4){c -(}
.                         estimates restore m1
.                         mchange_mi std_same_prop_UnEmpl "continuous" "mi"
.                         estimates restore m1
.                         mchange_mi std_same_prop_PhysProb "continuous" "mi"
.                 {c )-}               
.         {c )-}
. {c )-}
{txt}
{com}. 
. if ("`model'" == "logit"){c -(}
.         svy: logit Suic i.St `model_eq', or 
.         estimates store m1 
. 
.         if (`model_version' <= 4){c -(}
.                 mchange `margin_demographics', amount(all) delta(100) statistics(all) decimals(7)
.         {c )-}
. 
. {c )-}
{txt}
{com}. 
. if (`model_version' == 5 | `model_version' == 6) {c -(}
.         estimates restore m1
.         margin_interact std_same_prop_Sex Female 5 `model'
.         estimates restore m1
.         margin_interact std_same_prop_AgeGrp4 AgeGrp4 5 `model'
.         estimates restore m1
.         margin_interact std_same_prop_Race5 Race5 5 `model'
.         estimates restore m1
.         margin_interact std_same_prop_BornUSA BornUSA 5 `model'
.         estimates restore m1
.         margin_interact std_same_prop_MarStat5 MarStat5 5 `model'
. 
.         if (`model_version' == 6) {c -(}
.                 estimates restore m1
.                 margin_interact std_same_prop_UnEmpl UnEmpl 5 `model'
.                 estimates restore m1
.                 margin_interact std_same_prop_PhysProb PhysProb 5 `model'
.         {c )-}
. 
. {c )-}
{txt}
{com}. 
. 
. log close 
      {txt}name:  {res}<unnamed>
       {txt}log:  {res}/N/project/suicide_study/pnas_replication/results/log/mi_3.smcl
  {txt}log type:  {res}smcl
 {txt}closed on:  {res}20 Aug 2020, 16:47:48
{txt}{.-}
{smcl}
{txt}{sf}{ul off}